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Abstract 

It is shown how to construct many-particle quantum-mechanical spectra of 
particles obeying multispecies exclusion statistics, both in one and in two 
dimensions. These spectra are derived from the generalized exclusion principle 
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multiplicity formula. 
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I. INTRODUCTION 



Exclusion statistics was introduced by Haldane Q by postulating a "generalized exclusion 
principle" , according to which adding a particle of species b to a system reduces by g a b the 
number of single-particle states available for each subsequent particle of species a. Thus, 
the single-particle Hilbert space dimension is 

D a = D®- y E9ab(N b -5 ab ) , (1) 

b 

where is that in the absence of particles and N a is the number of particles of species a; 
the many-particle dimension is taken to be 

w = tx Cg. + N. - 1)1 

K l N.'.(D. - 1)! 1 > 

The quantities g a f, are called mutual statistics parameters and constitute the statistics matrix 
G = {g a b '■ a,b = 1, . . . , s}, where s is the total number of species. 

Starting from these formulas, it is possible to calculate the thermodynamic quantities: 
the distribution function (at least implicitly) and the equation of state (at least in the 
form of a virial expansion); see Refs. @,[| for the single-species case and Refs. || for the 
multispecies case. In particular, when the exponent a in the single-particle dispersion law 
e(p) = eo + Ap a is equal to the dimensionality of space, all the G dependence of the equation 
of state is contained in the second order of the virial expansion |^|. 

The question to be addressed here is the following. What is the many-particle spectrum 
of noninteracting particles obeying exclusion statistics in a potential with a given single- 
particle spectrum? Understanding this is an important step towards solving a more general 
problem of finding multispecies integrable models that realize exclusion statistics. Thus far, 
there are two known kinds of such models: (i) the single-species one-dimensional Calogero 
model (and the Sutherland model ||); (ii) the single-species P,|X0f and multispecies 



11, U lowest Landau level (LLL) anyon models. The spectrum of the Calogero model can 



be viewed either as that of interacting bosons/fermions or of free particles obeying exclusion 
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statistics [ 12], 13| , the 1/x 2 interaction being traded for a "renormalization" (fractionization) 



of quantum numbers, which can be equivalently regarded as implied by the generalized 
exclusion principle. (In some way, this is reminiscent of the equivalence of Aharonov-Bohm 



interaction and anyon statistics |14[].) The same structure of the spectrum occurs for the 
LLL anyon model, which can be also generalized to many species. 

In the present Letter, we will point out a simple and natural way of constructing mul- 
tispecies many-particle spectra (in a harmonic potential) for a symmetric statistics matrix 
{dab}, that are to be viewed as a realization of exclusion statistics in the sense that they 
have the generalized exclusion principle built in and lead to the thermodynamics implied 
by Eqs. (|l|) and (0). After reviewing the one-dimensional case, we make the construction in 
two dimensions, where the underlying idea is to introduce the fractional exclusion principle 
"along one degree of freedom" . 

II. ONE DIMENSION: CALOGERO MODEL 

We start with recalling the known single-species one-dimensional case. The single-particle 
energy spectrum in a harmonic potential with frequency u is 

E n = (n + ~) u , (3) 

where n = 0, 1, 2, . . . . An A-boson state is characterized by a sequence of quantum numbers 
{rij : j = 1, . . . , A} such that Uj+i > nj, and its energy is E{ nj y(0) = J2j=i E n .. The N- 
boson partition function is 

z " (o) ^S^ [ -^» (o)1= i^£w-ii - (4) 

where x = (3uo. The spectrum of the Calogero model, 



E {n 3 }{g) 



A/ 1\ A(A-l) 



w, (5) 

where g(g — 1) is the coupling constant of the inverse square interaction, can be obtained in 
terms of free particles obeying exclusion statistics with parameter g, in the following way: 



Regard n in Eq. (0) as a continuous variable and postulate that every particle k "pushes" 
every particle j residing higher than k, by g units up, in terms of n. This is understood as g 
states being excluded for particle j. Mathematically, this means that the quantum numbers 
are "renormalized" : rij — > hj , with 

JV 

™i = n i + ~ h k)9 , (6) 

k=l 

where 9 is the step function. With n,j+i > rij, the solution to this equation is hj = rij+(j—l)g 



12|1 ; that is, the lowest particle remains in its place, the second lowest one gets shifted by 
g units up, the third one by 2g up, etc. On the other hand, there is no interaction now, so 
the many-body energy is a sum 

TV 

£ {n . } (<7) = $> fl . . (7) 
i=i 

By summing (^|) over j, one sees that renormalizing the n/s increases the energy by gu per 
each pair of particles. Hence, Eq. (^) follows immediately. In particular, for g = 1 one obtains 
the fermionic spectrum. The partition function is Z N (g) = exp — |iV(iV — l)gx Zjv(O), and 
the thermodynamics is the same |3|-§| as obtained starting from Eqs. (|l|)-(0). 

III. A MULTISPECIES GENERALIZATION 



A state of s species of bosons is characterized by a set of quantum numbers {n a j : a = 
1, • • • , s, j = 1, . . . , iVj with n aJ+1 > n aj , its energy is E {naj} (0) = E«=i E^i E naj , and the 
partition function is 

Z Nl ,.. Ns (0) = f[Z Na (0)- (8) 

a=l 

We will assume the statistics matrix to be symmetric: g a b = gba — an important restric- 
tion, which will be discussed later. The appropriate generalization of (|^) is (cf. |15| ) 

s N b 

n a j = n aj + ~ n bk )g ab . (9) 

6=1 k=l 



4 



This means that a particle of species b "pushes" one of species a by g a b units up whenever 
the latter resides higher than the former in the final state. Obviously, in contrast to the 
single-species case, the order of n's may not coincide with that of n's, so that Eq. (|9|) has to 
be solved in a self-consistent manner. However, assuming again the free-particle expression 
for the energy 

8 N a 

%(G) = EE^. (io) 

a=l 3=1 

the form of Eq. implies that, just like in the single-species case, each pair of particles, 
one of species a and the other one of species b, contributes g^uj to the energy. Hence, 
r s V ' 1\ , °N a (N a -l] 



E{n a3 }{G) 



EEk + 9 9aa+ NaNtg* 

La=lj=l V Z/ o=l Z a,b=l 

a<b 



(11) 



This is essentially the same as the LLL spectrum of multispecies anyons (with charges 
of the same sign) [pTif4|. Moreover, a semiclassical treatment of that problem Q leads to 
Eq. @. The problem itself is two-dimensional, but the LLL restriction renders it effectively 
one-dimensional by reducing the phase space. In the genuine one-dimensional case, however, 
it is not known what Hamiltonian would be a multispecies generalization of the Calogero 
one and possess the spectrum ([□]). At least two possible such generalizations, one fT6 



involving and the other one [[17] not involving a three-particle interaction, were considered 
in the literature, but neither of those serves our purpose. 



IV. THERMODYNAMICS 

Be that as it may, Eq. (|ll]) by itself is enough to derive the thermodynamic functions. 
From the grand partition function H = ^2n 1 ...n s Zi 1 ' ' ' z ^ s Zn 1 ...n 8 [where z a = exp(/3[i a )\, one 
can deduce the cluster expansion and the dimensionless analog 

of the virial expansion In 2 = E^...^ A kl _ ks N^ 1 ■ ■ ■ N^ s /Z^'"^' 1 , using the relation N a = 
z a (dliaE/dz a )] the coefficients A kl _„ ks are just numbers. If there is a well-defined volume, 
the usual virial expansion (the equation of state) (3P = J2kx...k a fl fci...fc s Pi 1 • • ■ p k s s follows, by 
\nE = (3PV, N a = Pa V. 

5 



The partition function that one obtains from Eq. (|ll|) is 

Z Nl ... Ns {G) = exp 



• N a (N a - 1) • \ 

2^ 7) 9aa+ }^ N aN b g ab \X 

a=l 1 o,6=l / . 



J N!...N, 



(0) 



(12) 



a<b 



This yields the cluster coefficients in the thermodynamic limit, 



b kl ... ks {G) = F kl ... ks {G)x 



-i 



(13) 



F kl ... ks (G) = (-1) — - — ^ 11 ^/'V ///-/- : ( 14 ) 

lln=l pi, 9i, j=i 

P S _1,9 S _1=1 
Pn<Qn 

in the last sum, all the pairs (pj, qj) must be distinct, and each number from 1 to s must figure 
in the set {pi, . . . ,p s -i, qi, ■ ■ ■ , q s -i} at least once. All the numbers k%, . . . , k s are assumed 
to be different from zero (otherwise, an obvious invariance with respect to the renumbering 
of the species and the identity i 7 fc 1 .../ Cr0 ...o = F kl ___ kr should be used). For example: 

F i0 o(G) = ^(1 - 4<7 U )(1 - ^ n )(l - ^ n ) ; 
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i 71 32i(G) = i[l- (3gn + 2g 12 +g 13 ) 
o 



[1 - (35-12 + 2#22 + 923) 



X (3#12#13 + 2^12^23 + 913923) ■ 



To compare, the general formula for particles obeying exclusion statistics in a D- 
dimensional harmonic potential, implied by Eqs. ([!]) and (fj), is Q 



b kl ... ks (G) = (k 1 + --- + k s y u f kl ... ks (G) Z 1 , 



(15) 



and for a symmetric G it turns out that f kl ... ka (G) = (&! + ••• + k s )F kl ... kg (G). Since in the 
thermodynamic limit Z\ = 1/x, the system at hand does exhibit exclusion statistics. The 
dimensionless virial coefficients are 



A 2 = (20H - l)/4 , A u = g 12 ; 
A k = B k -i/k\ for k > 2 , 



(16) 



and all the others vanish. Thus, they only depend on the statistics at the second order 
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The system considered is, in fact, thermodynamically equivalent to a system of particles 
with linear dispersion e(p) = Ap in a box (because the single-particle spectra are equidistant 
in both cases, and the renormalization will alter them in the same way) of a volume related 
to the harmonic frequency in such a way that the single-particle partition functions Z% in 
both cases match: 1/x — V/nAf3, or V = ttA/u. The virial coefficients of the latter system 
are a kl ... ka = A k ,.. Je .{nA/3) kl+ - +k '- 1 . 



V. TWO DIMENSIONS 

The single-particle spectrum in a two-dimensional harmonic potential is 

E mn = {m + n + l)u , (17) 

where m, n = 0, 1, 2, . . . . The quantum numbers m and n correspond to the two degrees of 
freedom (say, x and y motion). The idea is to implement the fractional exclusion principle 
"along one degree of freedom", i.e., with respect to one quantum number. Thus, split up the 
spectrum into "sectors" , labeled by m, each sector being the spectrum of a one-dimensional 
oscillator with a constant (m + l/2)u added to the levels, and postulate exclusion statistics 
in such a way: Particles within one sector behave like they do in the above multispecies 
one-dimensional model, while particles in different sectors do not influence one another. 

Formally, it looks this way. A state of s species of bosons can be characterized by a 
set of pairs {(m a j,n a j) : a — 1, . . . , s, j — l,...,N a } such that either m a j+i > m a j or 
(m a j + i = m a j and n a ,j+\ > n a j). That is, the next particle of the same species is put either 
into a higher sector than the previous one or into the same sector but not lower. Define a 
set of numbers {N a i : a = 1, . . . , s, 1 = 1,..., / max }, where M a i is the number of particles of 
species a in the Zth from below nonempty sector and Z max is the total number of nonempty 
sectors. One has N a = ■Afai- Introduce now a quantity M. a i = J2\'=i-Afai'', then for any 
a, the Z-th nonempty sector contains those particles (aj) for which A4 a ,i-i + 1 ^ j < Mai- 
Eq. (g), in accordance with the aforesaid, is replaced by 
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s N b 

n a j = n aj + E E 9(™aj ~ ^bk)^m a] m bk 9ab 
6=1 k=l 

whereas m a j = m a j. The free-particle expression for the energy 

s N a 



(18) 



a=l j=l 



(19) 



now yields 



E {(m aj ,n aj )}(G) = E 



1=1 Lo=l 



E-^ U + £ + E E 



a=l j=M a ,, 



+ E o 9aa + E NalNugab 



a=l 



a,b=l 
a<b 
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(20) 



where mi is the absolute number of the Ith nonempty sector, so that m a j = mi with / such 
that M. a ,i-i + 1 < j < -Maz- The quantity under the outermost sum can be recognized as 
Eq. (|TT| ) applied within the l-th sector, plus a constant displacement determined by mi. For 
example, let in the case of Bose statistics there be particles of species 1 in the states (0,1), 
(2,4), (3,1); particles of species 2 in the states (0,0), (2,1), (2,2); particles of species 3 in the 
states (1,2), (2,2), (3,2). The energy of the corresponding state with exclusion statistics will 
be (3 + g 12 ) + (4) + (21 + 2g 12 + g 13 + g 22 + 2g 23 ) + (11 + g 13 ) = 39 + 3g 12 + 2g 13 + g 22 + 2g 23 . 
The partition function takes the form 

OO OO OO Imax ( / S \ 



z Ni ...nXg)= EE E •■ 

{Mai} "li=0 m2=mi+l mi 



E 



max 'max 



n n%a(o) 

-i+l 1=1 l \a=l , 



x exp 



E a/;, (mi + + Kl(Af ; 1 l) 9aa + e jw*g*)x 

.a=l ^ ^' 0=1 ^ a,6=l / . 

a<b 



(21) 



The outermost sum is over all possible distributions of particles over sectors, with l Tl 
determined by such a distribution. 

The cluster coefficients that one obtains from this are 



K...k s {G) = (fci + • • • + k s )~ l F kl ,„ ks {G)x 



(22) 



with F kl ...k s (G) as in Eq. (TTJ]). This is but again the general formula for D = 2, so that 
exclusion statistics is present. The virial coefficients do not exhibit any specific pattern. 



S 



It is possible, then, to obtain the cluster coefficients of the same system in a box of 
volume V rather than in a harmonic potential P,|IT|, |n||T9H : 

Vr 2 

bi... ks (G) = -^(h + ■■■ + k s )b kl ... ka (G) (23) 

(the dispersion law being quadratic, e{p) = p 2 /2m, and A 2 = 2nf3/m). These coincide, up 
to a common factor, with the ones from Eq. (0), and the dimensionless virial coefficients 
are again those given by Eq. fll6f) . (In fact, they will be the same for any system with 
exclusion statistics with a constant single-particle density of states.) The dimensional virial 
coefficients are a kl ... ks (G) = A kl ... ks (G)X 2 ^ + - +k ^ . 



VI. DISCUSSION AND CONCLUSIONS 

Essentially, what has been shown here is that (i) for a one-dimensional harmonic poten- 
tial, "excluding g ab (single-particle) states" as required by the generalized exclusion principle 
has to be understood as renormalizing the quantum numbers so that the sum of those for a 
pair of particle a and particle b increases by gab] (ii) in two dimensions, the same procedure 
has to be followed within effectively one-dimensional subsets of the spectrum, corresponding 
to one degree of freedom. In fact, we have also verified that the same thermodynamics is 
obtained when in the two-dimensional problem, the sectors are made up of states with the 
same angular momentum, which appears to be a better choice since it preserves rotational 
invariance. It is a plausible assumption that any division by sectors is actually good, and that 
the same will be true even for the three-dimensional problem [where Ei mn = (l + m + n + ^)uj 
and the sectors would be labeled by / and m]. 

It is, however, not clear at this point how to construct the spectrum if the matrix g a b is 
not symmetric. In deriving Eq. (|H|), o r (POP, it was implied that 8(x — y) + 6(y — x) = 1 
always. However, putting #(0) = 1/2 and using it for a nonsymmetric matrix does not lead to 
the same thermodynamics as obtained from Eqs. ([I]) and (|2[). Perhaps some generalization of 
Eq. (|9[), or (JT^), is necessary in this case. Apparently, there is a relevant physical example: 
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Interpreting numerically found few-electron spectra in terms of fractional quantum Hall 
effect quasiparticles and comparing the multiplicities to Haldane's formula [f20| seems to 
hint that the system of two species of anyons with opposite charges in the LLL should 
exhibit exclusion statistics with gx% = —gi\- To prove (or disprove) this directly by finding 
the spectrum, i.e., by generalizing Eq. (|TTD, remains an open problem. 

For a symmetric statistics matrix, both the spectrum and the thermodynamics of the 
conjectured integrable models realizing multispecies exclusion statistics are now known, 
which provides one with an essential information helping to search for such models. 
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